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Abstract 



flj ■ The separate study of quark and gluon jets is vital for the interpretation of 

,^ , multiple variables behaviour observed in both high-energy hadron and heavy- 

ion collisions in the present and future experiments. We propose a set of jet- 
energy dependent cuts to be used to distinguish between quark and gluon jets 
experimentally based on a Monte-Carlo study of their properties. Further, we 
p^ I introduce the possibility to calibrate these cuts via gamma-jet and multi-jet 

events, which represent clean production channels for quark and gluon jets, 
respectively. The calibration can happen on real data and thus, reduces the 
tlT' I dependence of the method performance on Monte-Carlo model predictions. 



1 Introduction 



Jets are objects produced in hard scatterings of colliding particles. Experimentally 
5^ I we can observe jets as showers of high momentum particles. The character of these 

showers is determined by the fragmentation properties of the original parton; quark 
or gluon. In QCD, quarks and gluons carry different color factors [1]. This factor is 
proportional to the probability of a parton to radiate a soft gluon. Gluons have more 
than twice the color factor as quarks and as such are expected to form broader and 
higher multiplicty jets with softer fragmentation function. 

Apart from these differences, the gluons are expected to contribute significantly 
to the baryon production as compared to quarks [2]. All these differences naturally 
must be demonstrated in the particle spectra observed in an experiment. Previous 
experiments with e~^e~ [3] and pp |1] collisions studying the properties of different 
parton types have qualitatively proven these expectations. 

From the heavy-ion perspective, the study of fragmentation properties of quarks 
and gluons becomes important for understanding unexpected observation from RHIC 



explained by different phenonienological models (e.g. coalescence [5], jet flavor con- 
version [6]), which incorporate the above mentioned differences. 

Our aim is to perform a systematic study of the baryon and meson production 
inside quark and gluon jets. For this, we need to identify the jets flrst. In the following 
we introduce a data driven method to distinguish quark and gluon jets and make the 
study of their properties experimentally feasible. 

2 Method description 
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Figure 1: Panel (a): < AR{90%) > as function of jet's energy for R = {0.4,0.7}. 
Lines represent the fits in separate energy intervals. Panel (b): Table with fit param- 
eters. Indicated are the fit functions for different energy intervals, where exp corre- 
sponds to AR{90%)caic = exp{A + Bx Ejet} and poll to AR{90%)caic = A + Bx Ej^f 

The experimental treatment of quarks and gluons in proton-proton collisions seems 
challenging since we only can observe the final state hadrons together with the un- 
derlying event. Therefore, our observation is restricted to the experimental definition 
of a jet in terms of jet finding algorithms. Further, the study of jet properties based 
on their identification may be biased by our prior expectations incorporated into 
Monte-Carlo models. 

However, experimental data offer the possibility to distinguish between quark and 
gluon in an unbiased way, by observing their properties in channels, were we are 
certain of the origin of the jet. Such channels are the multi-jet and gamma-jet events, 
sources to gluon and quark jets respectively. Observing the properties of jets in such 
events can help us to identify the leading jets in others. 

The study was performed on simulations done using the Pythia 6.4 Monte-Carlo 
generator with the settings of PerugiaO tune [7]. For testing 4 data samples were 
generated, each containing 1 milion events. The samples involved 3 sets of events 
with hard scatterings, divided based on the two leading jets into gluon-gluon (GG), 
quark-quark (QQ) and quark-gluon (QG) samples. In order to study quark jets, 
gamma-jet sample was created (7Q). To reconstruct the jets, we used anti-kx jet- 
finding algorithm [8] for 2 jet-size parameters; R = {0.4,0.7}. In order to design 
the cut, we chose to compare the subcone size, which contains 90% of jet's energy - 



Ai?(90%), after the tracks have been sorted in distance from the jet axis. The method 
introduced is performed in two steps. 

1^* step From each event we select two leading jets and measure their Ai?(90%). 
We plot it as a function o jet's energy (see Fig. [H panel (a)). We fit < AR{90%) > 
and so we obtain AR{90%)caic (see Fig. [1], panel (b)). 

2"*^ step Next we reconstruct jets in multi-jet and gamma-jet events in order to 
obtain samples of gluon and quark jets. From multi-jet events we select all but 
the two leading jets, from gamma-jet events we take the jet at 180° ± 30" degrees 
w.r.t. the gamma. We measure AR{90%) measured for the selected quarks and gluons 
and check how it is distributed around AR{90%)caic- We obtain a distribution of 
DR = AR{90%)caic — AR{90%) measured plottcd in Fig. [21 Based on this distribution 
we choose a DR cut to be applied to the leading jets (Table [I]). 
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Figure 2: DR variable distribution for R = {0.4, 0.7} and different jet energy intervals 
for MC quark and gluon jets compared to DR distribution of quark and gluon jets 
obtained from gamma-jet and multi-jet events respectively. 



3 Discussion 



We presented a possibility to distinguish between quark and gluon jets experimentally. 
The main advantage of this method lies in the fact, that the cut on the variable 
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Table 1: Cuts to select quark and gluon jets based on gamma-jet and multi-jet dis- 
tributions of DR in Fig. |2l 
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Figure 3: Performance plots. Efficiency and purity as function of Ejet- 

introduced can be calibrated on real data without the necessity to rely on Monte- 
Carlo information. This way, our selection is not biased by our prior expectations 
about the differences of quark and gluon jets. 

The performance of the method was determined by the efficiency and purity of 
the identification of the two leading jets in an event. As can be seen in Fig. El for 
quarks it achieves efficiency and purity up to 60% in the higher energy interval. For 
gluons, efficiency of the selection is better (constant, 60% for higher energy interval), 
although the purity drops slightly with energy. The rise of purity for quarks and the 
drop for gluons is however expected since quarks are to form harder jets than gluons. 
The method performs better for bigger jet size, namely R = 0.7 (see Fig. |3l panel(b)). 

The cuts were applied for two energy intervals. In Fig. [2] we see, that the gamma- 
jet and multi-jet DR distribution does not overlap well with the Monte-Carlo quarks 
and gluons in the lower energy interval, especially for R = 0.4, which distorts our 
selection and worsens the performance. We need to be aware of this and apply a more 
detailed treatment of this region. 



Acknowledgements 

I would like to thank my supervisor Peter Levai for cooperation and the support of 
my ideas. The presentation of this work was possible thanks to the funding provided 
by ELTE and the Hungarian OTKA 77816. 



References 

[1] R. K. Ellis, W. J. Stirling and B. R. Webber, Camb. Monogr. Part. Pliys. Nucl. 
Phys. Cosmol. 8 (1996) 1. 

[2] S. Albino, B. A. Kniehl and G. Kramer, Nucl. Phys. B 725 (2005) 181 
[arXiv:hep-ph/0502188] . 

[3] P. Abreu et al. [DELPHI Collaboration], Z. Phys. C 70 (1996) 179. 

[4] A. A. Affolder et al. [CDF Collaboration], Phys. Rev. Lett. 87 (2001) 211804. 

[5] V. Greco, C. M. Ko and P. Levai, Phys. Rev. Lett. 90 (2003) 202302 
[arXiv:nucl-th/0301093] . 

[6] W. Liu, C. M. Ko and B. W. Zhang, Int. J. Mod. Phys. E 16 (2007) 1930. 



[7] P. Z. Skands, Phys. Rev. D 82 (2010) 074018 |arXiv:1005.3457l [hep-ph]]. 

M. Cacciari, G. P. Salani and G. Soyez, JHEP 0804, 063 (2008) [ arXiv:0802.lT89l 
[hep-ph]]. 



